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ABSTRACT 


The  velocity  of  the  Interface  between  liquid  and  vapor  phases  when  a 
subtly  superheated  liquid  changes  to  vapor  is  found.  The  temperature 
distribution  in  the  liquid  ahead  of  the  evaporation  front  can  then  be  determined 
by  a  numerical  integration. 
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length  scale  factor  defined  in  equation  (20) 
slope,  dp/dt,  of  the  saturation  cuz*ve 
speed  of  sound  In  the  liquid 

position  (i.e.  x  co-ordinate)  of  interface  at  time  t 

velocity  of  interface 

thermal  conductivity 

latent  heat  of  e/aporation 

mass  rate  of  evaporation  per  unit  aurea 

pressure 

time 

temperature 
fluid  velocity 
space  co-ordinate 

x-g(t)  «  space  co-ordinate  measured  from  interface 
thermal  diffusivity 

y/A  »  dimensionless  distance  from  interface 

dimensionless  temperature  defined  in  equation  (9) 

density  of  liquid 

density  of  vapor 

o{  t/A  -  dimensionless  time 

dimensionless  velocity  of  interface  defined  in  equation  (14) 


SUBSCRIPTS 5 

o,  refers  to  conditions  at  interface,  liquid  side 
s,  refers  to  equilibrium  saturation  conditicm 
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INTRODUCTION 

This  study  is  an  attempt  to  solve  the  problem  of  the  motion  of  a 
plane  front  in  a  fluid  medium  across  which  a  change  of  phase  t^.kes  place. 
In  particular,  we  assume  that  an  infinitely  long  tube  is  filled  with  a 
liquid  at  a  uniform  pressure,  p(oo),  and  at  a  uniform  temperature, 

T(oo),  slightly  above  the  saturation  temperature  corresponding  to 
p(oo)«  The  liq\iid  is  thus  in  an  xinstable  state. 

We  assume  that  at  some  point  in  the  tube  (x^O)  some  vapor  is 
formed,  and  that  we  then  have  an  evaporation  front  moving  into  the  super¬ 
heated  liquid.  As  fluid  crosses  this  front  it  is  converted  into  a  vapor 
at  equilibrium  saturation  conditions.  The  energy  to  supply  the  latent 
heat  of  evaporation  is  assumed  to  come  by  conduction  through  the  liquid. 
Hence  the  liquid  will  be  cooler  at  the  interface  than  at  infinity. 

We  let  X  «  g(t)  be  the  location  of  the  front  and  g'(t)  its  velocity. 
The  determination  of  this  function,  g'(t),  is  the  principal  object  of  the 
investigation. 


ANALYSIS 


The  temperature  In  the  liquid  must  satisfy  the  one-dimensional 


diffusion  equation  in  a  moving  substance, 
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in  the  region  x2g{t),  ttO,  subject  to  the  boundary  condltims: 

=  -T  (2) 

(3) 


T(\^)  =T(oo) 

T(oo^t)  =  T(oo) 


(4) 


In  boundary  condition  (2),  which  expresses  the  fact  that  the 
latent  heat  of  evaporation  is  supplied  by  conduction  throuj^  the  liquid, 
the  functions  g(t)  and  m(t)  are  unknown.  If  we  assume  that  the  moving 
front  converts  the  liquid  into  stationary  vapor,  then 


from  which 


=  fv  5'^^)  =  f 
=  0  - 


(5) 

(6) 


At  212®F,  ■  7.5  X  10“^ j  hence  it  is  permissible  in  considering 


the  effect  of  motl<»i  of  the  liquid  to  take 

^0=  5  Vo 


(7) 


The  curved  boundary  x  g(t)  can  be  removed  by  choosing  a  co-ordinate 
measured  from  the  interface.  That  is,  letting  y  -  x-g(t),  equation 
(1)  becomes 


dt 


(8) 


The  tenpereture  gradient  will  be  large  only  near  the  boundary  y  ■■  0, 


and  here  Qi  -  g'(t^  is  small  because  of  (?)•  Hence  the  convective 
term  in  (8)  is  neglected., 

The  problem  is  reformulated  in  dimensionless  quantities  by 


writing 


T-T^eo) 


(9) 


where  T^  (oo)  is  the  equilibrium  saturation  temperature  corresponding 
to  p(oo).  We  then  have  the  standard  diffusion  equation 


a 


(10) 


(the  subscripts  denote  partial  differentiation)  to  be  solved  in  the 
region  Z’i  0,  with  the  boundary  condltims: 


©0;  0)  ~  0 

(11) 

e(oo,  r;  =  0 

(12) 

(0,  T)  .  ^(r)  (13) 

The  boundary  function  ^(Z*)  in  (13)*  still  to  be  determined,  is 
the  dimensimless  form  of  the  velocity  of  the  Interface,  g'(t),  as 
is  seen  from  (2),  (5),  and  (9). 


>4L  fv  d(t) 
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The  boundary  value  problem,  equations  (10),  (ll),  (12),  (13), 
is  most  conveniently  solved  by  ^>plying  the  Fourier  cosine  transform 
on  the  variable  7^ ,  Let  the  transform  be  denoted  by 

f  6(v,r)  df 


For  tho  boundarj  oondltions  of  tho  proiont  probloa  tho  p«rti«l  difforoixtial 


•quatlon  (10)  transforu  to 

~  <t>(r)-  \  © 

with  the  solution 
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It  can  b«  shown  that  the  iaTsrelon^of  (15)  Is 

~  f  ^  0Cs)('r-s)  ^  a/5 

giving  the  tenpsrature  distributicm  in  toms  of  the  motion  of  the 
intsrfaoe. 

Proa  (16)  we  find  the  tenpsrature  at  the  interface  to  be 

</>/5)  (17) 


Since  themodjmaaic  equilibriun  is  established  at  the  evaporation 
front,  the  boundary  tempsrature,  7^  (2^)  “  T(0,  ?'),  given  by  (17), 
will  determine  the  pressure  at  the  interface  p^(?’).  If  the  range 
of  temperature  variation  is  siaall  enou^,  we  can  take  the  relationship 
between  pressure  and  temperature  at  saturation  to  be  linear.  That  is. 


p^-p(oe)+b 

1 

We  now  assume  that  the  evaporation  front  moving  into  the  originally 
motionless  liquid  sends  a  wave  into  this  liquid  in  which  the  fluid  velocity 
and  the  pressure  are  related  as  in  an  acoustic  medium. 
p(x,t)  -  p(oo)  ♦  pcu(x,t) 

In  particular,  at  x  ■■  g(t) 

p^-.  p(.,)  r  fc  3'») 

%LA 


(19) 


If  «•  now  ohooM  our  length  aemlo  A  to  bo 

fcK  ,  . 

* "  r^Lb 

equating  (18)  end  (19)  providee  the  following  integral  equatim  for 

the  notion  of  the  interface: 

,T 


Equation  (21)  is  eolved^by  an  iterative  technique  in  idiieh 


(21) 


with  •  1«  the  Iteration  can  ^  carried  to  any  order  with  only 
elementary  integrals  being  encountered.  The  result  is  the  following 
power  series: 

'P-0*ZtiT‘i-i7U . . ) 

Aft«r  some  niAnipulatlon  the  series  can  be  reco^ized  ae  the  e3q>an8lon  of 
er/c  (22) 


idiere  erfc  denotes  the  complimentary  error  function  defined  by 

/✓ 


r*  I 

orfc  (  Z  )  “  1-  erf  (  ?  )  -  1- 1.  /  .  “ 

Wri  ^ 


c 

RESULTS 


Fig.  1  is  a  graph  of  equation  (22).  Fig.  2  shows  some  typical 
temperatiu*e  distributions  obtained  by  numerical  evaluation  of  the 
integral  in  equation  (l6). 


Fig.  I  Velocity  op  Interface  vs.  Tir-ig 
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